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abstract 

This article deals with a single period inventory model with inventory returns and special sales. The demand is 
assumed to occur in a uniform pattern during a planning period say, t p _ Both non- deteriorating items and 
deteriorating items are considered for the discussion. 
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I. INTRODUCTION 

The problem of returning or of selling the inventory excess to optimal stock level is considered for an 
order level inventory system, when the demand is probabilistic in nature. The optimal stock level of an 
inventory system is smaller than the amount on-hand, and then the question of selling excess inventory arises in 
many business organizations. That is in any wholesale or retail business the demand of particular product 
decreases due to launch of new product which is cheaper or superior. This is so happened in the context of 
Pentium V computers. That means the sales of Pentium V. This type of context can be materialized when a new 
budget suggests the price increase or due to any other market fluctuations the stock-on-hand of the product 
becomes more than the optimal stock level. In such instances, the optimum amount to be returned or sold, if any 
should be determined by balancing the losses due to various costs involved in the system. 

Naddor [1] has considered this problem for EOQ inventory system where as Dave [2] has extended for 
an order level inventory system. We now develop a single period inventory model with probabilistic demand for 
non- deteriorating and deteriorating items. Inventory models for deteriorating items have been considered by 
various researchers in the recent past, including Shah & jaiswal [3], Aggarwal [4], Nahamias [5 - 9] .In all these 
models the convention problem of ‘how much’ ‘and ‘when to stock ‘was discussed under varying assumptions 
of inventory systems. 

Choosing appropriate hypothetical parameter values one can carry out numerical illustration. The 
applicability of the model requires specific parametric quantification with different distributions. Even with 
simple functions the cost function becomes complex in nature.In analysing such inventory systems, it is always 
challenging to estimate the demand distribution. Most of the researchers use parametric approaches for 
Estimation of the demand distribution.Usually a gamma distribution cf.Moomen 1962[10], Burgin[ll]; Das[ 1 2] 
a lognormal distribution. (Beckmann and Bobkoski [13]; Tadikamalla[14] ) and other distributions. Since, with 
small data sets, it is difficult to reject a null hypothesis of any hypothesis of any of the above standard 
distribution; it is interesting to study a non-parametric approach. Strijhoch and Heuts [15] have used parametric 
and non-parametric approximation for the lead time demand distributions and demonstrated that the Kernal 
density (Non-Parametric approach) estimation approach has considerable advantage over the parametric density 
estimation in an (Q,s) inventory model with a cost criterion. A separate simulation study on these lines will be 
carried out highlighting the computational convenient over parametric approach. The pertinent results and 
analysis will be communicated elsewhere 

In this study probabilistic inventory model with inventory returns for non -deteriorating items when the demand 
will be a random variable X, during the period t p . 

II.MODEL ASSUMPTIONS AND NOTATIONS: 

We adopt the following notations in carrying the mathematical model. 

(i) X: Random variable representing the demand during the period with known probability density 

function f(x). 

(ii) The scheduling period is a prescribed constant say t p 

(iii) Replenishment rate is infinite, replenishment size is constant. 

Lead time is zero. The fixed lot size q p raises the inventory at the beginning of the scheduling 
period to the target level S. 
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(iv) Shortages if any are to be backlogged. 

(v) The system starts with an amount of Q units out of which only P units are to be retained, after returning 

or selling. The problem is to determine the optimal value say P" of P. 

(vi) The inventory holding cost Ci per unit time, the shortage cost C 2 per unit time and the returning or 

selling cost C 4 per unit and known and constant during the planning period t p . 

III.MODEL WITH NON-DETERIROTATING ITEMS: 

In this context the period T starts with initial inventory level of Q units and final inventory level of 
zero. Since (Q - P) units are returned (or) sold, the retained P units will be consumed by the time say ti (ti < T) 
due to demand. During the period (T - ti) the optimal order level system will operate. Two typical situations 
may arise in the system depending on the relative values of demand x during a cycle and the quantity retained P 
units. 

Case 1: x < P i.e., demand size less than or equal to the quantity retained. 

Let Qi (t) denotes inventory position at time t (0 < t < ti). Then the differential equation governing the system 
during t, is given by 

i Qi M =-x/t lJ o<t<t 1 (l] 

The solution of the above equation using the boundary conditions 
Qi (0) =P and Qi (ti) = 0 gives 

Qi (t) = [P-x/t (e‘-l)] e" 1 , 0 < t < * (2) 

To determine average number of units In(x) carrying inventory per unit time, we assume that Qi (t) is 
approximately linear in t, so that 

IilOO =;[(J L C0)+(J 1 Ct L ]] =; PCl+e-^-^Cl-e'O (3) 

Case 2: x > P i.e., P < X < Q x . Here we note that maximum quantity that can be returned (or) sold, if any is Q 
i.e. the optimal value of P must be less than or equal to Q. In this case, suppose that the system carries P units 
and these P units will be consumed during ti period and during ( ti ,T), the order level inventory system will 
operate. Then it can be seen that the inventory level Qi (t) over (0, T) is given by 

QiCtD = [{P - - 1)} e" 1 ] 0 < t < * (43 

= ^(*-0* <t <T 

Since Qi (0) = p and Qi (ti) = 0, equation (4) gives 

* = 108(1+3 (53 

Further since ti < T, we have 

x > Pti / (e* 1 — lj = P* (say) (6) 

the average number of units per unit time during the cycle are 

■»-£->* Ol p> 

Here we use the series form of exponential and logarithms in the above equations (4.3.6) and (4.3.7) we get 
Iu =£andP* = P (S) 

From equations (4.3.3) to (4.3.8) the expected average amount of inventory is 

K(P) = ^ [C 4 J^tP - x) f(x) dx -h C L I u (x) ttxjdx +■ C 4 (x - P) fix) dx -f (x)f(x)dx F 

(T-OkCS) ......... .(9) 

Where K(S) is the expected cost function as given by Naddor[l] . In order to reduce the complexity of the 
various expressions involved in the model we use series form of exponential terms and ignore the terms of 
second and higher order terms of P we get 

p 

l/ 1 ^ j"^P j" f j" ^ P^ \ 

K(P) = — [C 4 J (P - x) fGOdv -H CjJ -(2 -*3f60dx + J xftodxj + cJJ (1°3 

0 0 D 0 

The optimum P, = 0 yields 
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k(pO = c 4 



P P <3 

j" r FT r p CP CP 

+ c L J ffcOdx - c L J — f&Odx + c L J -fGO dx - TftP) 

0 D P 


= 0 - - (113 


The integrals involved in equation (10) and (11) can be evaluated explicitly for some specific probability density 
functions f(x) only; for those functions we can obtain optimum value of P (say) P° by solving (11) and then 
minimum total expected cost of the entire system can be obtained by substituting P = P° in equation (10). This 
is what illustrated in the following example. Let the demand density function is given by 

f(x3 = c[0.002e -<MMHx ] ,X > 0 


= 0 otherwise (12) 

Where C I 0.002 is used to normalize the above function as density function 
For this density function equation (11) becomes 

q l [q0,002e- <MNK< > - G.0G4Pe- M ® zP + - 2e" MI}zP + l] + j[l - e " MMP (0.002P + l)] + 

Ci p[ e -MMP e -MMQ] + £i£ C0.002) 2 e- OJHHP - ^ (0.002) 2 e-* 1 ™ 1 ’ = 0 . . (13) 


The solution of the above equation gives optimal quantity to be retained say P0 of P. However, explicit solution 
in P for the above function is not possible. Hence, to find the optimal value of P say P°, we use Newton- 
Raphson iterative method for solving the equation (13). 

IV. COMPUTATIONAL RESULTS: 

We have endeavored to determine what effect varying in the parameter C 4 i.e., disposable cost or 
special sales over the changes of optimum cost. Such changes can take place due to uncertainties in any business 
context or they might be influenced by the decisions maker himself. To examine these implications of the 
model, consider a hypothetical system with the following parameter values. Ci = 0.65; C 2 = 5.04; R = 2004; 
T = 0.5; and Q = 4800. All the parameters are expressed in consistent units per unit time. For different values 
of C 4 , we have determined the optimal quantity to be retained and associated costs are summarized in the 
following table 


C 4 

ti 

p o 

K(P°) 

S 

K(S) 

0.20 

0.104758 

272.7947 

4056.292 

500 

2882.321 

0.22 

0.115979 

306.1074 

5682.961 

500 

2882.321 

0.24 

0.126517 

338.1883 

7609.305 

500 

2882.321 

0.26 

0.136181 

368.2986 

9716.512 

500 

2882.321 

0.28 

0.144946 

396.1952 

11917.03 

500 

2882.321 

0.30 

0.152833 

421.7822 

14135.35 

500 

2882.321 

0.32 

0.159901 

445.1115 

16317.61 

500 

2882.321 

0.34 

0.166229 

466.3219 

18429.25 

500 

2882.321 

0.36 

0.171900 

485.594 

20450.61 

500 

2882.321 

0.38 

0.181582 

519.0852 

25473.23 

500 

2882.321 
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We note that, from equation (13) the optimum value of P° of P increases with an increase value of C 4 . 
This is obvious since the higher the retaining or selling cost, the smaller should be the quantity to be returned. 
Thus the optimum value of P is sensitive with respective to the parameter C 4 


V. LINEARITY ASSUMPTION AND ERROR ANALYSIS: 

In deriving the model, it has been tactically assumed that the depletion of inventory Level I (t) of the 
system is linear in t, over the cycle time. In fact, I (t) is a non-linear function of t and the linearity assumption 
over estimates the average amount in inventory, in this case average amount to be retained and consequently the 
total cost of the system. This will induce an error in the cost function. However, if this error is too small in 
magnitude then the assumption of linearity is justified. We now present an error analysis and show that incurred 
error is significant. Note that, if the linearity function is not made then the expression for average amount to be 
retained in inventory per unit time during the cycle is given by 


Ifi (t) = I ; 0 T In (0 dt = f Cl - e T ) - x (14) 

Li ti 

Hence the error due to linearity assumption is 

eI 1L (t L ) = 1^© - I^Ct) (15) 


Here 1, l(r) is given in (4.33) and the relative error in total cost REK (T) is given by the ratio 


racrft- 'l — K(ti)-K E (t i j _ CfCttj) 
EECt lJ " K E ttO " K E (ti) - 

where K(tj = C L Ij_ L (x) fix) dx 


( 16 ) 


K E (t L ) = C L ^>(-1 - e T ) - x]fO0dx 

- Ci[/ 0 P ' Wfix) dx - J. P x fix) dx] 

Here, K E (t L ) can be retained for the given model by replacing I 1L (t) with I* L (0 in equation (10). By 
taking = Re 0.56 per unit per unit time and ranging C 4 = 0.28 (0.02) 0.38; different values of K(t L ) (using 
equation 10), K E (T) are obtained. The corresponding relative error in total cost (REC) are obtained and 
summarized in the following table. 


Table 1: Error Analysis 


C 4 

P 

K(tJ 

K E (tJ 

REK 

0.28 

396.1952 

27.98553 

15.86942 

0.004228 

0.30 

421.7822 

32.72707 

18.7021 

0.004086 

0.32 

445.1115 

37.38957 

21.51662 

0.003989 

0.34 

466.3219 

41.90605 

24.26878 

0.003922 

0.36 

485.594 

46.23575 

26.92975 

0.003611 

0.38 

519.0852 

54.26194 

31.91789 

0.003443 
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From the above table, it can be observed that the magnitude of the relative error in total cost is very small which 
shows that the linearity assumption is justifiable. 

VI. MODEL WITH DETERIORATING ITEMS: 

The probabilistic model is developed with the same assumptions as those given in section 2 except 
that a constant fraction O(say) of the on-hand inventory deteriorates per unit time. Depending on the values of P 
and x of X, the following cases arise: 

Case (i): When x < P 

Let Qi (t) denote the inventory level of the system at time t (0 < t < t L ) when the demand x unit 

X 

occurs during the time ti then during that period the demand rate is — and hence differential equation governing 
the system is 

i Qift) + QiCO = =* , Co < t < t L ) (17) 

Ut ti 

The solution of (17) gives 

Q L (t) = [p — 7 "j~ fexp(0t) — l}] exp (— Ot) (0 < t < t L ) , 0 < 0 < 1 . (IS) 

Since x< P, we must have Qi (t L ) > 0 (or) x < P 1 ^/ {exp (Ot i - 1 ) } = P“(say) 

The number of unit, that deteriorate during the cycle is 

D,(x) = P-x-Q,(t,) = (l — e - * 1 ) (p + — x ,X < P ...(19) 

To determine average number of units I LI (x) carried in inventory per unit time we assume that Q L (t) is 
approximately linear in t, so that 

111 M = f [Qi CO) + Qi M = i P[1 + e- et i] -^[l- e -HX<P. (20) 

Case (ii): When x>P: 

In this case, suppose that the system that in inventory during (0, t L ) an during (t L ,T) the optimal 
order level system will operate; then it can be seen that the optimal inventory level Qi(t) over (0, T) is given by 
QiCtO = [ p - - l)] e _9tl ,0 < t < tj_ = [1,(5) + Ij (s)](t - t,Xt, < t < t, (21) 

Where i±(.s) and I 2 (s) can be had from Shah and Jaiswal [3] since Qi Ct L J = 0, (4.3.21) gives 

t, = Jlog(l (22) 

Further, since t t < T, we also have in this case 
x > OPT / {e 91 - 1} = P 

The no. of units that deteriorate during the cycle are 

D 2 60 = P-^log(l t^),x >P (23) 

The average amount in inventory per unit time is 
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=W = 5i l08 ( 1+! r)' ,:>f ’ 


(24) 


From (19) and (24), we have 

D(PJ = Jq D L (x}f(x)dx + D 2 (x)f(x) dx 

= Jr Pdt t (1 - — - x)f(x) dx - Jp — f(x)dx 


.(25) 


From (4.3.20), (4. 3.24) and (4.325), the expected total cost of the system is given by 
p fp Q Q 

k(P) = y j(P - x)f(x}dx 4- j 4- J I u (P}f(x)dx ^ j (x - P)f(x) dx 


c 

+ T 


f DWfOO + j D, Cx) /Cx) 


ax 


■f (T - O^CsO (26) 


Where K(S) is the expected cost in the order level system and is given by 
K(S) - (^0 4- c L j (s - -x^ /(x)dx 4- ~s 2 (.c& 4- c t + c 2 9T") J 

^ sfl-yST) 


/Of) 


dx 


+ ■ 


r Cx-s) ! s*&T r 

) 2x dx - T“ Cci + 2c ^ J 


s 3 9T 


/ Cx) 

~zr d * C 27 ) 


Optimum value of i ‘ of S, in the solution of the following equation 

(cB -\ci_9T + c L + c 2 )F [j(l - ; 9T)] + (c6 4- c 2 9T + c L 4- c 2 ) S ^_^rfdx - V0r( Cl 4- 


2c 2> dx-c 2 = 0 


.(28) 


For the derivations of K(S) in (27) and S in (28) see Shah & Jaiswal [3]. The optimum value P D of P is the 
solution of the following equation. 


Pi = 0 

=> K (P) = j 


f PT0 1 r P : l ! e : -PTs] f P 2 T9 (PT9 \) 

J ftxJdx + — 1\P) + Cl J [1 + ^ — |/(x3dx+|j>+ — (— -lJj/QO 

u Ld ) 

Q Q n 2 

Ch [TO f PT9\ c 4 f c p [ 2PTB / T&\ 

+ 7 l J T fVOdx -\ F -~) fn] + T l -l rtoW + TJ l ~ (i~t) 

P P 0 

Q 

T9 f T9\ PTB f 2 PTB 

~ —]f{x)dx 4- [PTB (l - — J - — ]/(P) - j —ftxydx + PTBf(P) (29) 


In practice, for obtaining the optimal value of P° of P, we are supposed to evaluate various integrals 

Sn ( jo 

of "'V ' =0. It is not possible to obtain explicit expression for some of the integrals. To overcome this 

difficulty we have used the expansion form of the exponentials and logarithms in the cost function (26) by 
ignoring the second and higher order levels of6 under the assumption thatO < T. The integrals involved in (29) 
can be evaluated only for some specific functions f(x) only, for such functions we can determine optimum 
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quantity to be retained say P" of P using equation (29) and the corresponding minimum cost can be obtained by 
substituting P° in (28). This is illustrated with the following example. Let the demand density function be 

/Or) = ~x 2 exp(— x) >x > 0 = 0 otherwise 

For this demand density, equation (29) becomes 

tftP) = h\(F') F K 2 (P') + K 2 (P) + A t 4 (P) +K 5 (.P) F K S (F) F (T - tjKCS) .. .. (30) 

Where 

K t (P) = |p[/>{-P*V p * - 2P x fi- pK 4 2(1 - e”P x 3} 4 P x * B - py 4 3P X * e" p * - fi{-P*e-P y 4)l 
+ e- p,t )}] .(313 

Px 1 ^' 


KJp) 


(, „„ , Px'e-n P\T8)\ . , 

= Cl p|(l- e - p j+ (e- p -l) 

p 2 -T 1 

4F*3 - l}-^ {e-^il + PO - 1} (323 


tfgOO = 7 [— Q 3 e pX — P x 8 pX — 3{— q 2 g p * \- P** e ** 4 2{— qs 4 F* 0 p * — 9 ^ + 9 p }] ... .. (333 

K 4 (P) = y [? 3 8-S - P^ 9~ P 1 - 3Q 2 8~* - 3F^ - 9-** 4 6{~Q 8~* 4 P* + (*“«- 0 "^)} 


- F + F* 2 e~ pI + 2 {{-q9~< 4 P x 9~ pI ) 4 (- 0 ^ 4 0- p 3}}] (343 




P l T9 _ u x\ t 
2 


-e- p (P* 4 13 4- (PTe3 ! (e"P x - 13 4 
F 6 (F3 =^[e-ni4q3 - e-P x (l + F*)] ... 


4 ' 2 


4 i P : e -P* 4 . ^ 


-(353 


(363 


AT (53 = ( C tf - ^ Cl OT 4 c L + cj) (s - 


4 - [4 ted 4 c L +cj3 (S 4 33 - c l S9T(S 4 63 - 2cjOT(2S - 3)] exp(-S3 (373 

£> 

and the optimum order level of S s ofs is a solution of 

— - c l 0T F c t F c 2 j — -Cc 4 eOT F e L F e 3 )(S F 2) — CjiGTCS 2 F4S — 6) — 2c- 3 0S(2S — 5)]exp (— S) 

= 0.. . .....(30) 

For derivation of K(S) and SI see Shah and Jaiswal [3]. The optimum quantity to be retained P° of P is the 
solution of 

K'tPj^) F K'(P 3 ) F K'(P a ) F K"(P 4 ) = 0.. (39) 


K ( Pl 3 =^[{-P x2 a- pK - 2PV* + 2(- e P s + l)} + £ j^V^} 

+ ^[ q3 e- q - Pe-P 1 - 2{e- pK (l + P x 3 - (q 413^}]. 


.(403 


K'(Pj3 = y [{-P x e" py - 2P* e p * 4 2 (l - e py )} 4 3 P 2 (OT3 3 -PTe - 2 P* e" py + 2(l-e" py ) 

4 [ P(Te3 J - FTe}P x e-P 3 ‘] (413 

K'(P 3 3 = ^ [7 {-q 3 - e- q 4 P 1 e“ pr 4 2 { e - py Cl 4 P*) - (q 4 l3}} - ( P “ ?)] '" f42 3 
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k(pJ = [Preti - (1 + ^ (i - e~n -2 { - 

p^g-P* -sp^e-P* -|-2(l - e - plt ) + ip x3 e - plt [™_^] 

PT9{Cq -I- IV - e -P * (l 4- P*)} 4- 7 PTe(P x V plt ;)] (43] 


VII. DISCUSSION 

In this article we have attempted to derive an analytical solution to the optimum quantity to be 
retained in the context of probabilistic demand. Both non-deteriorating and deteriorating cases were considered. 
However, a closed form solution could not be obtained. To handle this type of situation one can use a non- 
parametric approach in estimating demand distribution. In all these cases the Kernel density estimation approach 
proposed by Strijbosch and Heuts [15], can be used. This suggestion will enhance the further scope of this 
research problem. 
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